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We study nonequilibrium evolution in a self-interacting quantum field theory invariant under 
. . . space translation only by using a canonical approach based on the recently developed Liouville-von 

CO ' Neumann formalism. The method is first used to obtain the correlation functions both in and beyond 

' the Hartree approximation, for the quantum mechanical analog of the model. The technique 

, involves representing the Hamiltonian in a Fock basis of annihilation and creation operators. By 

separating it into a solvable Gaussian part involving quadratic terms and a perturbation of quartic 
}^^' terms, it is possible to find the improved vacuum state to any desired order. The correlation functions 

, for the field theory are then investigated in the Hartree approximation and those beyond the Hartree 

■ approximation are obtained by finding the improved vacuum state corrected up to O(A^). These 

' correlation functions take into account next-to-leading and next-to-next-to-leading order effects in 

0^ , the coupling constant. We also use the Heisenberg formalism to obtain the time evolution equations 

for the equal-time, connected correlation functions beyond the leading order. These equations are 
derived by including the connected 4-point functions in the hierarchy. The resulting coupled set 
\l \ of equations form a part of infinite hierarchy of coupled equations relating the various connected 

^ . n-point functions. The connection with other approaches based on the path integral formalism 

is established and the physical implications of the set of equations are discussed with particular 
emphasis on thermalization. 
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PACS numbers; ll.lO.-z, ll.lO.Ef, ll.lO.Wx, ll.15.Bt 
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p i' I. INTRODUCTION 

^ ^ In the past few years, a lot of attention has been focused on the investigation of classical and quantum fields 
■ evolving out-of-equilibrium. Such interest is physically well motivated because the very early history of the universe 
^ ■ provides many scenarios where nonequilibrium effects may have played an important role. The reheating of the 
• ^ , universe after inflation , formation and growth of domains in any generic spontaneous symmetry breaking phase 
^ ' transition , the formation of topological defects 0, 0, IE IE Hi and possible formation of Quark-Gluon Plasma 
during the deconfinement transition or Disoriented Chiral Condensates during the chiral phase transition 9] are just 
some instances where the proper understanding of the physical process may crucially depend on our understanding of 
nonequilibrium quantum fields. The experimental accessibility of some of these phenomena, like formation of quark 
gluon plasma, made possible through heavy-ion colliders at RHIC and LHC, has also been a strong motivating factor 
behind the revival of interest in nonequilibrium evolution of quantum fields. The development of new theoretical 
techniques and the availability of more efficient computational resources have also made it possible to investigate in 
some detail the nonlinear effects which play such a crucial role in the evolution dynamics. 

Another important issue in this context is that of thermalization in closed quantum systems. Is it possible for 
macroscopic irreversible behavior to manifest itself starting from microscopic reversible (unitary) quantum dynamics 
[ift llll Il^ ? In other words, would it be possible for a closed quantum system to thermalize when it is perturbed 
from its initial thermal state? Can the process of thermalization be adequately described within the mean field 
description? If not, what is the role of interactions (non-linearities) in bringing about thermalization? Is it possible 
to develop a consistent theoretical framework to address these important issues? Some of these questions have only 



'Electronic address: senguptaOphys.ualberta.caJ 
^Electronic address: khannaCSphys.ualbertaxaJ 
Electronic address: sangkim@kunsan.ac.kn 



2 



rece ntly begun to be addressed using newly developed theoretical tools for dealing with nonequilibrium quantum fields 

[laiiiiiiMiiiia. 

Till recently, the issue of thermalization of a closed quantum system was typically addressed by a separation of 
the system into a sub-system made up of the non-thermal soft modes (with longer thermalization time-scales) and 
an environment consisting of the thermal hard modes (which thcrmalize on much shorter time scales compared to 
the soft modes) 0, 0, . The subsequent interaction between the soft modes and the environment (often treated 
stochastically) leads to the eventual thermalization of the sub-system made up of soft modes. Mean field theory 
(Hartree approximation) has also been extensively and successfully used to study the dynamics of nonequilibrium 
quantum fields 0, 0, 0, 0| and has yielded valuable insights into the early-time dynamical behavior. However, 
since the mean field (Hartree) approximation is essentially a linear approximation, it fails to correctly capture long- 
time dynamical behavior where nonlinear effects play a dominant role. Moreover, in the mean field approximation, the 
different field modes interact with a spatially homogeneous mean field with equal strength and therefore the effects of 
direct scattering (which is responsible for redistribution and eventual equi par tition of energy among different modes) 
are neglected. Hence thermalization of the system cannot be achieved [Tfil |. Recently, there have been some 
attempts ^B, '2l| to explore the possibility of thermalization within the mean field scheme by considering spatially 
inhomogeneous mean fields. The motivation was to see if the scattering of the fluctuation modes with the modes of the 
inhomogeneous mean field aids in thermalization. Salle, Smit and Vink poj proposed a density matrix describing an 
ensemble of pure Gaussian initial states. By defining the mean field to be an ensemble average over a set of Gaussian 
density matrices, and considering the interaction between the quantum fluctuation modes and the inhomogeneous 
mean field, they showed approximate thermalization is observed over intermediate time scales in the sense of particle 
distribution of low momentum modes approaching a Bose- Einstein form for a spontaneously broken theory. For long 
times, the particle distributions were found to tend towards a classical Boltzmann form. However, in a contrasting 
study, Bettencourt et.al. |2l| showed that by considering the dynamics in the presence of a spatially inhomogeneous 
mean field, but without ensemble averaging over a certain class of initial conditions, they found that the Hartree 
approximation fails to establish a thermal Bose-Einstein particle distribution even at late times. In a more recent 
paper Salle and Smit pointed out that the energy densities used in the simulations of [IJ were not large enough to 
ensure even approximate thermalization in the BE sense over the time scales observed. They found that by choosing a 
large enough energy density and a different set of initial conditions, BE behaviour was exhibited by the low momentum 
modes at intermediate time scales, even for a symmetric theory. However, even the approximate thermalization time 
scale was much larger in comparison to that of a spontaneously broken theory. In none of these works, quantum 
thermalization was observed for all modes and all energy densities, over the time scales observed. 

In order to understand the complicated process of thermalization of closed quantum field systems, one needs to go 
beyond the mean field (Hartree or leading order Large - N) approximation schemes. This has been attempted recently 
by many different groups 13 ll5ll23L l23| . More recently, time- reversal invariant equations for correlation functions 
of nonequilibrium scalar fields have been derived based on a 3- loop expansion of the 2PI effective action |2^ |2^ . 
This method has also been extended to 0(N) symmetric field theories by carrying out a systematic 1/N expansion of 
the 2PI effective action. The 3-loop approximation (for the single real scalar field case) and the next-to-leading-order 
(NLO) expansion for the 0(N) model incorporates the effect of direct scattering and shows promising evidence of late- 
time thermalization. A detailed analysis of the nonequilibrium dynamics shows the evidence of three distinct temporal 
regimes characterized by an early-time exponential damping, an intermediate hydrodynamic regime and a late- 
time regime manifest through exponential approach to thermal equilibrium. An alternative approximation scheme has 
been developed by Cooper and collaborators focusing on resummation methods based on the exact Schwinger-Dyson 
equations This scheme called the bare vertex approximation scheme (BVA) (first studied by Kraichnan |23) 

involves obtaining the exact Schwinger-Dyson equations and then neglecting the vertex corrections. All the methods 
discussed above are based on path-integral techniques within the Keldysh-Schwinger closed time path formalism 
and the basic dynamical equations can be quite complicated to derive. 

Our main aim in this paper is to make use of canonical formalism to study the nonequilibrium evolution of quantum 
fields. The Liouville-von Neumann (LvN) formalism was first applied (within the Hartree approximation scheme) to 
study the early-time growth of domains during a quenched second order phase transition [29i | . The LvN formalism 
which solves directly the quantum LvN equation, is another quantum picture [sol Isij besides the Schrodinger and 
Heisenberg pictures. Further, it is shown that the LvN formalism provides a convenient and powerful method for 
analysing nonequilibrium systems such as time-dependent oscillators ^30. .31. .32] and quenched phase transitions 
[2^ I33L Im . l35l |. We use this formalism to obtain the Gaussian vacuum and thermal evolution equations in the 
Hartree approximation. We develop a method for obtaining the improved vaccum state and explicitly show how to 
obtain equations for the first and second order coefficients which lead to the improved vacuum state correct to O(A^). 
We then derive expressions for the 2-point functions by taking the expectation value of products of field operators 
with respect to this improved vacuum state. This method clearly indicates that non-Gaussian effects first make their 
appearance at 0{X^). 
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An alternative canonical approach based on the Heisenberg formalism is then employed to obtain a set of nonequi- 
librium evolution equations for the correlation functions. After clarifying the relation between the LvN formalism and 
the Heisenberg-picture, we take the vacuum expectation values of the Heisenberg equations for all possible combina- 
tions of products of field operators and obtain a hierarchy of coupled equations for the ordinary n-point correlators. 
To obtain the dynamical equations for the connected equal-time correlators beyond the leading order, we make use 
of the method of cluster expansion which allows us to express the ordinary n-point correlators in terms of their con- 
nected counterparts. This method provides an alternative non-perturbative approach for going beyond the mean field 
approximation. We compare our results with those obtained in the literature using other methods 2 1. Mli . In this 
context we consider two recent approaches used to obtain the evolution equations for the partition function (3^ and 
correlators in a quantum mechanical model discussed by Ryzhov and Yaffe (RY)(2J|. We show that the canonical 
approach used in this paper yields the appropriate evolution equations, and thereby establish a connection between 
various approaches for obtaining the nonequilibrium evolution equations for equal-time correlation functions. 

The paper is organized as follows. In the next section we briefly outline the Liouville-von Neumann (LvN) formalism 
which will be subsequently used to study the nonequilibrium dynamics of the correlation functions in and beyond the 
Hartree approximation. In Sec. Ill, we describe in detail, a quantum-mechanical model of the anharmonic oscillators 
and use the LvN formalism to obtain the evolution equations for the vacuum and thermal correlation functions in 
the Hartree approximation. In Sec. IV, we develop the LvN formalism to study the non-Gaussian dynamics of 
the quantume mechanical anharmonic oscillator model. This section sets the stage for the application of the LvN 
formalism to study non-Gaussian dynamics in the more complicated field theory model. The issue of stability of 
the LvN method, in the context of the anharmonic oscillator, is discussed in Sec. IV. B. In Sec. V, the evolution 
equations for a nonlinear self-interacting scalar field theory are first derived in the Hartree approximation, both for 
the nonequilibrium as well as for the thermal equilibrium case. Nonequilibrium evolution beyond leading order is 
discussed in detail Sec. VI which contains the most important results of this paper. The LvN formalism is used 
to investigate nonequilibrium dynamics beyond the leading order in Sec. VI. A. In Sec. VLB. we make use of the 
Heisenberg formalism to obtain a hierarchy of nonequilibrium evolution equations for the connected n-point functions. 
A comparison between our approach and other methods used in the literature is carried out in Sec. VII. Sec. VIII 
contains a summary of our main results and a discussion of the physical implications of the results in the context 
of issues like thermalization in quantum field theory and phase transitions. The method of cluster expansion [s^ is 
outlined in Appendix A, where we also write down the set of equations for the connected correlators in configuration 
space. The equivalence between the RY method and the LvN formalism is established in Appendix B. 

II. LIOUVILLE-VON NEUMANN PICTURE FOR TIME-DEPENDENT SYSTEMS 

The systems under study in this paper have Hamiltonians whose coupling constants (parameters) depend on time 
explicitly. These systems describe nonequilibrium processes in the sense that the Hamiltonians do not give the correct 
density operators e^'^^^'^/Zjj. To properly find the Hilbert space and density operators for such time-dependent 
systems, we have to clarify the picture of the quantum evolution. Let the Hamiltonian be defined in terms of the 
Schrodinger operators for a quantum mechanical system 

Hit)^H{ps,qs,t), (1) 

and for a quantum field theory 

Hit)^H{nsJs,t). (2) 

Here the systems depend on time only through time-dependent coupling constants (parameters). In addition to the 
well-known Schrodinger, Heisenberg and interaction pictures there is another picture for such nonequilibrium systems 

MM 

Firstly, in the Schrodinger picture, the (functional) Schrodinger equation (in unit of /I = 1) 

Himit)) (3) 



u{tms (4) 



H(t)U(t). (5) 



^^l*(t)) 
has the exact quantum state 

determined by the unitary evolution operator 
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For the time-dependent case in contrast with the time-independent case, it is not easy to find the evolution operator, 
which is formally defined as 

U{t)^Texp J H{t')dt' , (6) 

where T denotes the standard time-ordered operator. Secondly, the Heisenberg operators 

OhW = UHt)dsU{t) (7) 

satisfy the Heisenberg equation of motion 

i^^duit) + [HH{t),dn{t)]=0, (8) 

where Hnit) is the time-dependent Hamiltonian Heisenberg operator. The Heisenberg operator Huit) is simply given 
by replacing the Schrodinger operators ps and qs by pu(t) and qii(t) in Eq. and -its and 4>s by nait) and (fiait) 
in Eq. (|2Jl. However, the explicit form of pait), Quit) and Tr-ait), 0H(i) requires either the exact knowledge of the 
evolution operator ^ in advance or the solution to Eq. Thirdly, we may introduce the Liouville operators 

di^{t)^u{t)dsuHt)- (9) 

Thus the Liouville operators evolve the Schrodinger operators backward in time . It follows then that the Liouville 
operators satisfy the quantum LvN equation 

i^6L(i) + [0L(i),ff(t)] =0. (10) 
We can show that any eigenstate of the Liouville operator 

ddmx,t) = x\^^,t), (11) 

where A is the eigenvalue of the corresponding Schrodinger operator Os, satisfies the Schrodinger equation (PJ. In 
fact, it follows that js^l 

|^(t)) ^ ^CAe*-^''*<*^'*l(**-^(*»'*^'*>|*A,i) (12) 

A 

where Ca's are constants. 

The essential idea of the LvN method Isol l3lL Is^ l33j | is that the quantum LvN equation provides all the 
quantum and statistical information of nonequilibrium systems. Technically the linearity of the LvN equation allows 
any functional of an operator satisfying the LvN equation to be another operator. Thus we may use some suitable 
operator Oh{t) to define the density operator p{t) = e~^^^^*^ / Zq for the time-dependent system provided Oi,{t) 
satisfies the LvN equation. In this sense the LvN method unifies quantum statistical mechanics with quantum 
mechanics. The LvN method treats the time-dependent, nonequilibrium system exactly in the same way as the time- 
independent, equilibrium one. Moreover, the LvN method can be applied to nonequilibrium fermion systems with a 
minimal modification 1391. 



III. ANHARMONIC OSCILLATOR IN THE HARTREE APPROXIMATION 

As a precursor to the investigation of (?5>'*-field theory, we apply the LvN method to a simple quantum mechanical 
model of anharmonic oscillators 

H = \p^±'^3:^ + ^S:\ (13) 

and derive the evolution equations for the coherent state expectation value of position and momentum variables as 
well as the subtracted 2-point correlators. The anharmonic oscillator with the lower sign is a quantum mechanical 
analog for the second order phase transition. All the time-dependent operators in Sees. HI, IV and V will denote 
Liouville operators, whose subscript L will be dropped. 
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The main idea behind the LvN method is to require the pair of invariant operators [32, |40j defined as 

a{t) ^i[u*{t)p-u*{t)x], a\t) ^ -i[u{t)p - u{t)x], (14) 

to satisfy the LvN equation i.e. 

^^ + [a,H{t)]=0, ^^ + [a\H{t)]^0. (15) 

Here u{t) and u*{t) are auxihary variables in terms of which the 2-point correlators will be expressed. The invariant 
operators may be made the annihilation and creation operators satisfying the standard commutation relation at equal 
times 

[a{t),a\t)]^l, (16) 

which leads to the Wronskian condition 

u*{t)u{t)~u*{t)u{t)^i. (17) 

All the other commutation relation vanish 

[a(i),a(t)] = [at(t),at(t)] = 0. (18) 

Equation (|14|) can be inverted to express the position and momentum operators x and p in terms of the annihilation 
and creation operators as 

x{t) = u{t)a{t) +u*{t)aHt), 

p{t) = u(t)a{t) +u*{t)a\t). (19) 
The coherent state is defined either as the eigenstate of a{t) 

a{t)\a,t) = a\a,t), (20) 
with a complex eigenvalue a or as the displaced state of the vacuum state given by 

°° n 

\a,t) = D^ia)\0,t) = e-^'^^^'T^M, (21) 

„=o 

where D is the displacement operator 

i)(a) =e-""'(*)+"*"(*l (22) 

The coherent state can also be found using the variational principle [4l| . 
The coherent state then leads to the expectation value of x, p, and p^: 

X = {a,t\x\a,t) = au(t) + a* u* (t) , 
p = {a,t\p\a,t) = au(t) + a* ii* (t) , 
{a,t\x'^\a,t) = x'^ + u* {t)u{t) , 

{a,t\f\a,t) = + u*{t)u{t). (23) 

The subtracted 2-point correlators defined below are then given by 

gxx{t) = {x^) -x^ ^ u*{t)u{t), 
9pp{t) = {f) -P^ = u*{t)u{t), 
9xp{t) = (xp) - xp = u*{t)u{t), 

9px{t) = g:^{t)^u*{t)u{t), (24) 

from which we obtain the evolution equations for the 2-point correlators 

9xx{t) = 9xp{t) + 9px{t), 
gpp{t) = u*{t)u{t) + u*{t)u{t), 

gxp{t) = u*{t)u{t) + u*{t)u{t). (25) 
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A. Correlation Functions in Coherent State 

The expectation value of the Hamihonian with respect to the coherent state (|20() leads to an effective Hamiltonian 

H^^ ^ {HU = ± + + \9^(±^ + + ^^9iT + 0{X'). (26) 

The auxiliary field u{t), after differentiating Eq. (|15|l with respect to x and taking the coherent state expectation 
value, satisfies the equation 

il{t)±uj^u{t) + {a,t\-^\a,t)u{t) = 0, (27) 

so does its complex conjugate u*{t). Inserting the expectation values obtained in Eq. (|23|) into Eq. (|27|l . one obtains 
the equation for the complex u(t) 

U{t) + + ^;g2 _^ ^(^^ ^ Q (^28) 

Using the above equation and the definition, Eq. (|24|l . of the subtracted 2-point correlators, we get the following 
coupled set of evolution equations for the subtracted 2-point correlators 

9xx[t) = gxp{t) + gpx{t), 

9ppit) = -(±w2 + ^x2)(5,j, + gp,)+0(A2), 

9xp{t) = gpp- g,,(±Lj^ + ^x^^ +0{\^). (29) 

It is important to note that, these equations are correct only up to 0{X^) except for the first one which is exact. 
This is due to the quartic term appearing in the potential, unlike the case of the simple harmonic oscillator where the 
quadratic (Gaussian) form of the Hamiltonian leads to the vanishing of all correlators greater than second order. In 
this case, even if one starts from a Gaussian state peaked at a; = for which all the correlators greater than second 
order vanished, the subsequent evolution of the coupled set of equations would induce the appearance of non- vanishing 
values for higher order correlators at 0{X^) and higher 24|- This will be discussed later when we obtain the evolution 
equations for a (/)^-field theory beyond the leading order. 



B. Correlation Functions in Thermal State 



For the case of the anharmonic oscillator in an initial thermal equilibrium with the positive sign for the unbroken 
symmetry, the evolution equations for the auxiliary variable and the 2-point correlators up to 0{X^) can be obtained 
by taking the expectation value with respect to the coherent-thermal state 

PCT = e^p[~P{naHt)a{t) + Sa\t) + S*a{t) + eo}], (30) 

ZcT 

where eo = ^^/2 + |(5p/(ri). As a{t) and a,''{t) satisfy approximately the quantum LvN equation. Eg. 1)15(1 . so does the 
density operator H3U|I . We choose SI to satisfy the gap equation for the unbroken symmetry 

n^ = „» + A. ,31, 

As discussed in |2^, the displacement operator (|22() with a = 6/ (ft) unitarily transforms the coherent-thermal density 
matrix to a thermal one 

D^ia)pcTDia) = ^ e,^p[-pna\t)a{t)] = pT. (32) 

By making use of the unitary transformation 

D\a)a{t)D{a) ^ a{t) - a, t)\a)a\t)t)(a) ^ a) {t) ~ a* , (33) 
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it is easy to show that 



(x^)cT = + coth 



where now 



(P'>CT = + COth(^)7r(i)M(t), 



X = (a, t\x\a, t) = —{au + a*u*), 
p — (a, t\p\a, t) = —{ail + q*u*). 



(34) 



(35) 



The expectation value of Eq. (|15|l with respect to the coherent-thermal state (|3(J|I leads to the equation for the 
complex u{t): 



Hit) 



±tJ^ + -X 



A 
2 



2 ' , -coth(^V*^^ 



2 / 



u{t) = 0. 



(36) 



The corresponding thermal 2-point correlators are then given by 



gTxx{t) = {x^)cT-X^ ^COth(^!^^U*{t)u{t), 
9Tpp{t) = (p2)cT-p' = C0th(^)r(i)M(t), 

gT2;p(i) = [xp) CT - xp = cofhl^^^ u*{t)u{t), 
g-Tpxit) = g^^pit) ^ coth(^^^u*{t)u{t), 
from which we obtain the evolution equations for the thermal 2-point correlators 

9Txx{t) — gTxp{t) + gTpx{t), 

grppit) = 



(37) 



u u 



gTxp{t) 



±.^ + ^.^ + ^coth(^) 



5Tpp — gTx 



(gTxp + gTpx) -f ©(A^), 
0(a2). 



(38) 



In the T = limit, Eqs. 1)37(1 and ((38|l reduce to the evolution equations ((28(1 and ((29(1 . Expressions for the subtracted 
2-point correlators can then be obtained by solving Eqs. 1(36(1 for u{t) perturbatively in powers of A. 

IV. ANHARMONIC OSCILLATOR BEYOND THE HARTREE APPROXIMATION 



In this section we apply the LvN method to the anharmonic oscillator, in Sec. Ill, to go beyond the Hartree 
approximation. Our stratagem is to represent the full Hamiltonian p3(l in the Fock space basis ((14(1 and to follow 
the standard perturbation theory by taking the quadratic terms as an unperturbed part and the quartic terms as a 
perturbation. For the convenience of computation, we express the representation in the normal ordering where all 
stand to the left of a. 

We divide H into a Gaussian part Hq, the quadratic part, and a perturbation Hamiltonian Hp, the remaining 
quartic part. 



H ^ Hr 



XHv 



where 



Hp 



1 

2 ' 

1 

4! 



±- 



: r : +• 



6A 
4! 



x^)g ■■ x^ ■■ +Eg, 



(39) 

(40) 
(41) 
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Here Eq is the vacuum expectation value 



2\2 
G- 



Using the normal ordered operators 

: {ua + u''a 



^ A:!(n-fc)! 



y*(n-fc)yfc^t('^ 



we obtain the expectation value with respect to the Gaussian vacuum state that is annihilated by a(t) 



1 r 

2 



• * • I z * 

U U±LU U U - 



A 



{u*u) 



(42) 



(43) 



(44) 



A few comments are in order. First, the separation of the Hamiltonian into the quadratic and quartic parts in Eqs. 
(|40|l and H41|) is reminiscent of the Caswell-Killingbeck method 47] , which separates the Hamiltonian into a solvable 
part and a perturbation. In fact, as we shall show below, the quadratic part is solvable via the LvN method even 
for explicitly time-dependent systems. Second, the quadratic part (|4U|I involves a term proportional to the coupling 
constant A, which makes any perturbation theory based on it reliable even in the strong coupling limit of A. This 
term is the same as the Hartree approximation, — > 6{q^)q^. As will be shown below, the wave function(al)s of the 
Hamiltonian H4U|) are the same as those from the Gaussian effective potential method by Chang and Stevenson [isf . 
The vacuum state is also the same as the Hartree approximation. The equivalence of the vacuum state between the 
LvN and Gaussian effective potential will be shown below (see also Ref. [29|'). 

The Hartree approximation in Sec. HI is equivalent to simply using the truncated quadratic part, Hg, as the 
unperturbed part. Separating the quadratic part Hq is the essence of the Hartree approximation H or the variational 
Gaussian approximation 48] . In fact, the perturbation XHp, though Hq contains a term of the same order as the 
perturbation itself. Then the truncated quantum LvN equation 



.da 



leads to the mean field equation 



[a, He 



Hit) 



0, 



9at 



_l_ 2 , * 
±UJ + —u u 



[at,7?G] =0, 



i(i) = 0. 



(45) 



(46) 



The mean field equation above can also be obtained by minimizing Eq in Eq. H44I) . which proceeds by varying with 
respect to u* , using Sii* /Su* = d/dt, and treating u* and u independently. The equal-time commutation relation now 
is guaranteed by the Wronskian condition of Eq. H17(l . The Gaussian vacuum state is annihilated by a 



a{t)\0,t)G = 

and the excited number states are obtained by applying 

at"(t) 



\n,t)c 



0,t)( 



These are the exact quantum states of the time-dependent Schrodinger equation only for Hq 



i^^\n,t)G = HQit)\n,t)G. 



(47) 



(48) 



(49) 



By fixing the time-dependent phase factor [32 
wave functions are given by 



and including the factor from c- number term in Eq. (|40|) , the harmonic 



«'G,„(a;) = 



1 



1 



1/4 



(2n+l)/2 



V22n! \ 2'Ku*u 



^/2^ 



exp 



(50) 



where Hn {x) are the Hermite polynomials and the last c- number term in the exponent comes from the corresponding 
c-number term in Eq. (|40|l . From now on we denote the wave functions (I50|l by \n, t)Q without any loss of generality. 
These states form an orthonormal basis of the Fock space: 



^,{n,t\m,t)Q = (5„ 



(51) 
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A. Beyond the Hartree Approximation 



To go beyond the Hartree approximation, we need to include the perturbation, which is now given by 

1 



(52) 



There was an attempt in Ref. 46i| to solve Eq. (|15|l for the full Hamiltonian including the perturbation H52(l by 
improving a and . Here we find the improved quantum states by directly solving the Schrodinger equation for the 
full Hamiltonian. The perturbation excites and de-excites any Gaussian number state \n, t)^. As {|n, t)Q\ constitutes 
a Fock basis, we therefore expand the exact quantum states as 



oo oo 



^'*)=EE^'C'«„(t)|m,i)G, 



where the lowest order coefficient is 



C(o) _ X 



(53) 



(54) 



Using the fact that any state, |m,t)G individually satisfies the Schrodinger equation H49|l for Hq, the Schrodinger 
equation ^ for the full Hamiltonian H13|l leads to the set of equations 



1=0 m=0 /=0 m=0 

Comparing the powers of A, we finally obtain a hierarchy of dynamical equations for the coefficients 

oo 



(55) 



(56) 



Another expression for Eq. H53|l may be obtained in a compact form using operators and a. For each fixed m, 
we first sum over I 



C„;™(t)=£A'CW„(i), 



(57) 



and then write any state, \m,t)Q as either an excited or a de-excited state of the given lowest order state \n,t)Q, 
which is realized by applying the creation or annihilation operators a certain number of times. Hence, Eq. H53|) can 
be written by introducing an operator, Ui as 



\n,t) = J2 C„;™(i)k,t>G = Ui[a\t),a{t)-t,X]\n,t)G, 



(58) 



m—0 



and the Schrodinger equation leads to 



i—Uiit,\) + [Uiit,X),HG]-\Hp 
ot 



Using Eq. H45I) and technically assuming that all time derivatives act only on c-numbers but not on the operators 
and a, we obtain an interaction picture-like equation for the operator U : 



i^Ui{t,X) = XHpUi{t,X). 
ot 



(59) 



We then obtain the formal solution 



\n,t)^Uiit,X)\n,t)G, (n = 0,1,2,. ••), (60) 
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where 



Ui(t,X) = Texp 



iX / Hpdt 



(61) 



Here T denotes a time-ordering for the integral and a^{t) and a{t) are treated as if they are constant operators. 

We now find the improved vacuum state up to any order either by solving Eq. H56|l or by acting with the operator 
in Eq. H61|l on the Gaussian vacuum state. For instance, the improved vacuum state to order is given by 



\0,t)[2] = |0,t)G + A^ cW„(<)|m,i)G + A2 ^ C^%(t)\m,t)o, 



m— 



m—0 



where the only nonvanishing coefficients are 



u*\t'), 



(62) 



(63) 



and 



-if 



4! 
3 



u*^t')u{t') f u*\t") 



-if 



1 



3%/2 
1 
4! 



t ft 



u\t') f u*^{t"). 



(64) 



The non-Gaussian nature of the vacuum state (|62|l can be exploited by calculating the kurtosis (higher moments). 
The two-point and four-point correlators with respect to the Gaussian vacuum state (|50|l are 



G(0,<|x^|n,i)G = 
G(0,t|£4|o,<)G = i{u*uf, 

whereas those with respect to the improved vacuum state (I62II are given by 

i(2)_2 , ^(2)*^ *2\ , /^(2)* , ^(2) 



[2](0,t|i;^|0,i)[2] = u*u + A2[V2(Ci;>2 + C)^^*u*^) + (C^J + + ^Cl^l* Cl^>^)u* u] + 0{X') 



(65) 



(66) 



+ 6V2{C^'^l{u*u)u^ + C^^l*{u*u)u*^) + (123C^^i*C^^i + SC^^^* + 3C^'o^)(m*u)2] + ©(A^). (67) 



Stability of the LvN Method 



The stability of the LvN method should be checked since it is a time-dependent perturbation theory. That is, any 
secular coefficient in Eq. (|53|l should be removed systematically to insure the physically meaningful solution. For 
that purpose, we compare the LvN method with the standard perturbation theory for the well-known anharmonic 
oscillator with unbroken symmetry (positive sign in Eq. H13|l ). In that case we find the solution to the auxiliary 
mean-field equation 



where 17 is given by the gap equation 



uit) 



1 



iflt 



2n 



UJ 



_X_ 
4fi' 



(68) 



(69) 



11 



Then the time-dependent wave function (|50|) is given by 

A 



*G,n(a;,i) = exp 



exp 



^i\^{n+ 1/2) - 
-ii n{n+ 1/2)- 



3202 
32^2 



1 fnV^' f ^\ 

— Hn\ Vilq \ exp 



V22nT\ TT 

X "^gAx) 



(70) 



where ^G,n(a;) denotes the harmonic osciUator wave function. 

Substituting the solution (|68|l into Eqs. (|63|) and 164() . we find the improved vacuum state corrected to 0{X^ 



|0,t)r 



exp 



A 



\ 2 32rj2, 
A V3A2 



1 + i 



A2 



25^6173 29 721^6 



|4)( 



29 • 31)5 
X2 

^ I 



A2 



27 . 3\/2176 
70A2 



|2)( 



212 .31)6' 



0(A3 



(71) 



(2) 

Note that Cq.q originating from the four quanta creation and the subsequent annihilation leads to a secular term 
increasing as t. This is not a drawback of the LvN method but just a consequence of the time-dependent perturbation 
theory searching time-dependent states. As the coefficient of |0,t)G is the first two terms of exp[iA2/(29 • 3ri5)f], it 
can be approximately absorbed into the overall time-dependent factor to C'(A2) 



exp 



A 



A2 



321)2 29 • 31)5 



(72) 



This factor coincides with the time-dependent solution to the full Schrodinger equation obtained from time- 
independent perturbation method (s^l, where the corrected energy at C'(A2) is 



n 

E[2] = ^- 



X 



A^ 



321)2 29 • 31)5 



(73) 



The higher order terms that come from the creation of even number of quanta and its subsequent annihilation of 
equal quanta or vise versa also contain secular terms proportional to powers of t depending on the number of such 
processes. All these terms will provide the correct energy to the Schrodinger equation. 

Another way to understand this phase factor and thereby secular terms is to use the formal solution 1)61(1 . The 
operator has an exponential form 



Ui{t, X) — exp 



-iX I Hp{t') + {-iX)^[ dt'H{t'), dt"Hp{t")]+0{X^) 



A c-number term from the commutator 

A2 rt 

exp 



(4!)2 



dt'u^it') / dt"u*Ht")[a\a^^] 





A2 


exp 


*29.3D5* 



(74) 



(75) 



is nothing but the phase factor H72|) . Now the time-dependent vacuum state to C'(A2) does not involve any secular 
term as shown 



|0,t)r 



A2 



|0)( 



A2 



27 . 3\/2D' 



-|2>( 



A 



V3A2 

25^61)3 29\/2D6 



|4)c 



27 . 321)6 ' 



70A2 

212 . 3^6 



|8)c 



0(A3). 



(76) 



We thus have shown that the seemingly secular behavior can be removed systematically by taking the proper time- 
dependent phase factor for the wave function. This phase factor yields the correct energy for the anharmonic oscillator 
(|13|l to any desired order. The idea of removing the secular terms of higher order corrections by absorbing them into 
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the corrected energy is equivalent to removing the secular terms by renormalizing frequency in the multiple-scale 
perturbation theory [s^ • We note that the LvN method proves very accurate because the lowest order vacuum state 
is a Gaussian state that extremizes the Hamiltonian and the corrected vacuum state is expanded in the Fock basis 
[so]. Further, the LvN method is a powerful tool to finding the quantum states for nonequilibrium systems that are 
undergoing phase transitions p9l |. For the field theory case, the elimination of secular terms can be achieved using 
multiple-scale perturbation theory and this aspect has been discussed in detail in [s^ . 

There is another kind of instability originating from the dynamics of the system itself. Phase transitions from 
an explicit symmetry breaking in time provide such a dynamical instability. The anharmonic oscillator H13|l whose 
quadratic potential changes signs from positive to negative can be a quantum mechanical analog for phase transition. 
After the symmetry is broken, the mean field equation 



2 I * 

-bJ + —U U 



u = (77) 



may have a period when the quadratic term (— i^-'^) dominates over the quartic one (Am*u/2). Then the u grows 
exponentially as u k e^^ j^phl] until the quartic term grows and becomes comparable to the quadratic one. During 
this period all higher order corrections C^'^ beyond the Hartree approximation grow exponentially as powers of u and 
u* . This dynamical instability ceases when the state reaches the true vacuum state and oscillates over it. Thus the 
dynamical instability for a limited period does not cause any serious secular behavior as for the static system. 



/)*-FIELD THEORY IN THE HARTREE APPROXIMATION 



Now we will apply the LvN method to the (/)^-field theory, but first work out correlation functions within the Hartree 
approximation in this section. The 0*-field theory to be considered in this paper has the Hamiltonian in D-space 
dimensions 



In 1 , -,n m -o A 



(78) 



where tt(x) — (j){x) is the conjugate momentum operator. We divide the quantum field (f) into a classical background 
and quantum fluctuations over this background, (f>{x) = 4>c{x) -\- 4>f{x)^ where the classical background (mean) field 
(j)c{x) is in general considered to be spatially homogeneous. Then the Hamiltonian can be decomposed as 



H{t) = H,{t) + Hf{t) + H;nt{t) + SHintit) 



(79) 



where 



Hcit) 

ijint(i) 
SHintit) 



4! 



2 , 1 ^V7;l ^2 I 22 I U 



2^^/ + 2(^<^/) +— n 



4! 



(80) 



Here He and Hi are purely classical and quantum parts. We have divided the interaction into Ifint with even 
powers of (/)/ and SHi^t with odd powers of ttj and 0/. The ipcix) is nonvanishing only for the case when symmetry 
is spontaneously broken. We treat individually and collectively the modes in the momentum space of quantum 
fluctuations and the inhomogeneous background field. The momentum modes are given by the Fourier transform of 
a field (operator) and its inverse transform 



F{x) 



[dk]F{k)e 



ik-x 



■ik-x 



(81) 
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where F denotes either (j)f{x) and Ttf{x) or (j)c{x) and 7rc(a;), and 

m = (82) 

Then the quadratic integral takes the form 

j d^xF^ix) = j [dk]F{k)F{-k). (83) 

The Fourier modes of 0/ and tt/ will be denoted by 

= y d''x4>f{x)e-"'-^, 
Uk = I ci^x7r/(x)e-*'=-^. (84) 



The hcrmiticity oi (f)f and nj implies that = and n[ = Il-k- 
The commutation relation of the fields 

[hix)^^fiy)] = ^H^-y), (85) 

leads to those of modes in the momentum space 

[h,,tik]=i{2Trf6{k' + k). (86) 
In terms of the annihilation and creation operators satisfying the equal time commutation relation 

[ak'{t),di{t)] = {27r)''6ik'-k), (87) 
the momentum space operators and Ilfc may be expressed as 

= fk{t)dk{t) + <fe(i)aLfe(i)> 
tlk = Mt)ak{t) + cp*_k{t)a''_kit)- (88) 
Here it is assumed that (p-k{t) = ipk{t) and 

^lit)^kit) - ^imkit) = i- (89) 
Then the annihilation and creation operators are also expressed as 

ak{t) = tivimk - ^imk], 

al{t) = -i[ip^k{t)fl-k-^-k{t)^-k]- (90) 

Note that the momentum space operators 4^ and Ilk are regarded as time- independent ones whereas ak{t) and a^(t) 
as time-dependent Liouville ones in the LvN picture. The Gaussian vacuum is the state annihilated by all dk{t): 

ak{t)\0,t)G=0, (91) 

or the product of the Gaussian vacuum state for each (t) 

\0,t)G = l[\Ok,t)G. (92) 

k 

In the Hartree approximation we consider only those quadratic terms from Hf and Hint whose Gaussian vacuum 
expectation values do not vanish. Then the interaction term has the Fourier modes 

^int = ^ / Idki] J [dk2] J [dks] J [dfc4](2^)^5(fci + fc2 + fcs + k4)Mki)Mk2)Kh. 

^ ^ J[dkl] J[dk2]{2wfd{kl+k2)Mkl)Mk2) J[dk3]K^-ks 

= ^<l>lix) I mh^.k. (93) 
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The quartic term in the Hamiltonian for the fluctuation field (j)f can be approximated as 
Then the resulting quadratic part takes the form 



He 



[dk] 



(94) 



(95) 



Under the field redefinition Ii±k = (life ± n-fe)/2 and $±fe — ($fe ± $_fc)/2, this Hamiltonian is equivalent to that of 
harmonic oscillators 



where 



(96) 



(97) 



We may identify the classical background (jjc and tTc either with the vacuum expectation values of (f) and tt or with 
the coherent state expectation values of ipf and iff, respectively: 



(98) 



= (./)(x))vac = {<Pf{x))cs = J [dfc]($fe)ese''=•^ 

TTcix) = (7r(a;))vac = (7r/(a;))cs = j [dk]{tlk{t))cse'^'"' , 
where (■ ■ • )cs denotes the expectation value taken with respect to the coherent state: 



(99) 



Now the equations for the auxiliary field variables and Lp*j^ can be easily obtained in the Hartree approximation 
by making use of the LvN equations for dk{t) and d\{t) for the Hamiltonian H96(l . Then the LvN equations 

^■^^^^')+[4(t),i7G(t)] = 0, 



dt 
. ddkjt) 
' dt 



+ [ak{t),HG{t)] = 0, 



(100) 



lead to the equations 



A 



A 



[dk']ipl,(pk' 



(101) 



2"" 2 

The equation for the classical background 4>c{t) is obtained from the effective classical Hamiltonian, Hc{t) + (i/int(i))G, 



0c(x, i)-V2,/),(x,t) + 



0c(x,i) = 0. 



(102) 



It is more advantageous to work with the evolution equations for correlation functions, rather than the field equa- 
tions. Firstly, the 2-point correlation functions in thermal equilibrium are related to the Bose-Einstein distribution 
function (for a quantum theory) and the temperature of the system (for a classical theory) in a fairly simple way; 
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thereby allowing us to use them as benchmarks to track the evolution of the system towards thermal equilibrium. 
Secondly, it is more convenient to make systematic improvements to the mean field description, by working with equa- 
tions for the correlation functions. As we shall see in Section VI, the mean field equations for the 2-point correlation 
functions are part of an infinite hierarchy of evolution equations for the connected, equal-time, n-point correlators 
and are obtained by truncating this hierarchy at the level of the 2-point functions. 
We define the 2-point subtracted correlation functions for the fields as 

gii{x,x';t) = {[4)f{x) - (j)c{x)][(j)f{x') - (l)c{x')])cs= {^f{x)$f{x'))^ 

g22{x,x';t) = ([7r/(x) - 7rc(a;)][7r/(a;') - 7rc(a;')])cs = (7r/(a;)7r/(x'))v 

gi2{x,x';t) EE {[$f{x) ~ (f>c{x)][nf{x') ~ nc{x')])cs= {(i)f{x)nf{x))^ 

g2l{x,x'-t) EE {[T^f{x)-TTc{x)][4>f{x')-(t>c{x')])cs={7^f{x)4>f{x')). (103) 

The subscript 'v' implies expectation value with respect to the vacuum state. Using the above expressions we get 



311(2; 


x';t) = 






522(2; 


x';t) = 






912(2; 


x'-t) = 


j [rffc]^^(t)^fe(t)e^'=-(- 


-x') 


521(2; 


x'-t) = 


j [dk]ipl{t)^u{t)e^^-^''^ 


^x') 



(104) 

From Eq. (|104|l . the equal-time 2-point correlation functions Gij{k,t) in momentum space can be defined through 
the Fourier transform 

g,j{x,x';t) = y"[dfc]Gy(fc,t)e*("-"'\ (105) 

where i, j 1,2. Taking the time-derivative of the 2-point correlation functions yields the following evolution equations 

gii{x,x';t) = gi2{x,x';t) + g2i{x,x';t) 



g22{x,x';t) 



[dk][rk{t)Mt) + v>imk{t)]e 



ik-{x — x ) 



gi2{x.x'-t) = [dkmit)Mt) + rMMW''-^^-'''^ 



(106) 



Substituting the expressions for (pk and (fl in Eq. I|106l) and making use of Eq. H105|l results in the following 
evolution equations for the correlation functions in the momentum space in the Hartree approximation up to 0{\^): 

Gii(fc,t) = Gi2(fc,t) + G2i(fc,i), 



G22ikA) 



A 



A 



+ ttC(<) + o3ii(0,t) (Gi2(fc,t) + G2i(fc,t)) + O(A^), 



2 , '^j,2/- ^ 



Gi2(fc,i) = G22(fc,t)- [^^. + ^0^(t) + ^gii(O,i) Gn(fc,i)+0(A2) 

one can define a quantity 7~^(fc) in tern 
which is found to be conserved in the Hartree approximation 



^ 2^... , (107) 

As shown in p^l23||. one can define a quantity 7~^(fc) in terms of the momentum space 2-point correlation functions, 



7-2(fc) = Gn{k,t)G22{k,t)~^[Gi2{k,t)+G2iik,t)f 



1 



(108) 



as a direct consequence of the Wronskian condition, Eq. H17|) . Issues of renormalization of this model have been 
investigated in detail in Ref. [i^. 

It is also possible to define thermal 2-point functions by taking the expectation value of the field operators with 
respect to an initial Gaussian thermal state 



1 



{^{x,mx',t))T - Tiie-f'«"^(j>f{x)(bfix')] 



(109) 
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where /?o is the initial inverse temperature of the system and 

oc 

Z EE Tr[e-'3«^<^] = ^K,t|e- 



/3of2fc(o+afc + i) 



nk,t). 



(110) 



The corresponding equations for the auxihary field variables, are again obtained from the LvN equations using the 
Hamiltonian (j96|l but with the expectation values taken with respect to the initial thermal Gaussian ensemble. This 
leads to the set of equations 



+ J^4>l{x, t) + ^\ I [dk]{2nk' + lM,ipk' 



Vlit) = 0, 
,{x,t) = 0, 



where nk{T) is the Bose-Einstein distribution for a field theory system at a temperature T ~ 1/Pq: 

1 



nk{T) 



(111) 



(112) 



Note that at T = 0, the above set of equations reduces to Eqs. (|101|) and (102). The expressions for the thermal 
2-point equal-time correlation functions are also easily obtained 



9ii{x,x';t) EE {(j){x)lj){x'))T 
ff^2(a;,a;';^) = (7r(a;)7r(a;'))T 
5^2(2;, a;'; t) = ((/)(a;)7r(a;'))T 
g^^{x,x']t) = {■k{x)^{x'))T 



[dk] coth 
[dk] coth 
[dk] coth 



^fc \ * ik-(x-x') 



[dk] coth 



2 



(113) 



from which it follows that 



Gfi(fc,i) = coth 

G22{k,t) = coth 

Gj2{k,t) = coth 

G2i{k,t) = coth 



[-^j^kV'k, 
I P^^k\ . 



ff3o 



Ufk'Pk, 



^k^k- 



We can define another conserved quantity for the thermal 2-point correlation functions 

h^r'ik) - G{,ik,t)GUk,t) ~ \[G^2ik,t)+GUk,t)]^ = J 



(114) 



(115) 



As in the case of the nonequilibrium evolution, it is straightforward to obtain the evolution equation for the thermal 
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2-point correlation functions in the Hartree approximation. Using Eqs. (|114|l and (|lll|l . we get 



Gi2(fc,t) = GUk,t) 



'^l + W + ^glii^, t)] Gfi(fc, + o{\^). 



(116) 



The above set of equations together with Eq. (|lll|l describe self-consistently the evolution of the thermal 2-point 
correlation functions in the Hartree approximation. It describes the evolution of the system initially in thermal 
equilibrium at a temperature Tq, after interactions are turned on at time t = Q. They reduce to Eqs. (|101|l and H1U7|) 
in the T = limit. It is also possible to obtain the set of equations Eqs. (|lll|l - (|116f) by taking the expectation value 
of the LvN equations, Eq. H100|l , with respect to the coherent-thermal state with Hq being the Gaussian (quadratic) 
part of the Hamiltonian obtained after Hartree factorization of Eq. H78I) , (but without decomposing the field into a 
classical background and fluctuations) using the same techniques employed in the example of the anharmonic oscillator 
in Sec. HI. The initial thermal ensemble can be described adequately by Hq and the solution to Eq. Hlll|) can be 
obtained as 



1 



-iClkt 



(117) 



where is given by the gap equation 



I [dfc'lT^coth 



(118) 



For t = 0, we then have (pk{0) = l/VSHfe. The initial thermal 2-point correlation functions are then given by 



Gfi(fc,0) = 
Gj^{k,0) = 



1 
2 



2 



1 



2 



(119) 
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It is useful to compare our results with those obtained in 12,11 144| , in the context of thermalization in classical field 
theories. After neglecting A-dependent terms, the solution of Eq. (|lll|l is (pk{t) = e"*"''*/ ^/2uJk- With this form of 
ifk, one can easily show that in the classical limit (h — > 0) the 2-point correlation functions reduce to 



Gj,{k,0) 


To 








= To, 


G{,{k,0) 


= Gl*ik,0) 



(120) 



iTo 



These useful relations, which arc valid for a classical thermal system with unbroken symmetry, have been used as a 
benchmark to study thermalization in classical field theory [23, |4J| . It is possible to carry out numerically, a spectral 
analysis of 2-point, equal-time, field and momenta correlation functions and track their dynamical evolution. The 
fiattening out (momentum independence) of the Fourier transform of the tttt equal-time correlation function would 
be an indicator of thermalization in classical field theory. 



VI. NONEQUILIBRIUM EVOLUTION BEYOND THE LEADING ORDER 



In this section we use two different approaches to discuss the nonequilibrium evolution beyond the leading order 
Hartree approximation. First, we use the LvN formalism to obtain expressions for the 2-point and 4-point functions 
correct to ©(A^) by including all quartic terms of the Hamiltonian and solving the time-dependent Schrodinger 
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equation as in Sec. IV. This method unifies both the LvN formahsm and the Schrodinger picture because the LvN 
formahsm is only used at the Hartree approximation and all the non-Gaussian contributions from quartic terms are 
found in the Schrodinger picture. The second stage is similar to the interaction picture, though all states are expressed 
in terms of the time-dependent Hartree basis. It is found that the 2-point functions have non-Gaussian contribution 
at order of 0{X^), confirming the result that the non-Gaussian effects appear only at 0{X^) Secondly, we make 
use of the Heisenberg formalism to obtain a set of evolution equations for the connected correlation functions beyond 
the leading order. We find that the set of equations form an infinite hierarchy, akin to the BBGKY hierarchy in 
statistical mechanics [45l | and incorporating next-to-leading order effects amounts to truncating the hierarchy at the 
level of the 4-point connected correlation functions. 

The key differences between the two formalisms lie in the fact that while the LvN formalism (like the Schrodinger 
formalism) deals with evolution of states (with the field and conjugate momentum operators being time-independent), 
the Heisenberg formalism is associated with evolution of field and conjugate momentum operators (with the quantum 
states being time-independent). However, even in the LvN formalism, certain operators like the annihilation and 
creation operators and any function of them satisfy the LvN time-evolution equation distinguished from the Heisenberg 
equation through a difference in sign. In fact, the LvN equations can be thought of as the backward time-evolution 
of the Heisenberg equations (38] . Furthermore, the LvN approach provides a perturbative method for going beyond 
the Gaussian approximation and is useful for calculating non-Gaussian effects on domain growth in theories with 
perturbatively small self-coupling constant '45|. On the other hand, the Heisenberg formalism provides a systematic, 
non-perturbative approach for going beyond the leading order Hartree approximation in studying the noncquilibrium 
evolution of quantum fields. Inclusion of connected n-point functions in the hierarchy of evolution equations is 
effectively equivalent to a loop expansion in powers of p^ . 



In the LvN formalism discussed in the previous sections, the time-dependent vacuum state is found approximately 
at the lowest order as the state that is annihilated by the annihilation operators of all momenta in the Hartree 
approximation which satisfy the LvN equation. Going beyond the leading order amounts to determining the proper 
vacuum state of the full non-linear, interacting field theory, which can be expressed in terms of a complete set 
of number states of the Gaussian Hamiltonian at the Hartree approximation. The Gaussian vacuum state can be 
improved to any desired order by including the perturbation part. The n-point functions are then obtained by taking 
the expectation value of the appropriate product of field and momentum operators with respect to the improved 
vacuum state. 

The normal ordered Hamiltonian is decomposed into the quadratic Gaussian part and the quartic perturbation 



A. LvN Formalism: Evolution beyond the leading order 



H — Hq + X''Hp, 



(121) 



where the Gaussian part is 



Hg 




(122) 



and the perturbation 



Hp 



^ 4! y id,ki][dk2][dk3][dk4]S{ki + k2 + ks + k4) (pki(Pk2Vk3'Pk4akiak2ak3ak^ 
+ 4:(p*_k,'Pk2'Pk3'Pkia^^kiO'k2ak3aki + 6(plkif-k2'Pk3Vk:ia'i^^a'i^^ak-,ak4 




(123) 
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We define the improved vacuum state as 

= E E E ^'C'«^,...K'"2,..;i)G (124) 

i=0 ni,n2,--- fci,fc2,--- 

where we have used a concise notation to represent ni particles with momentum ki, 712 particles with momentum k2 
and so on, for the state and the coefficient. The summation over momentum is for modes which have non- vanishing 
particle numbers. The subscript 'G' refers to the Gaussian vacuum state. For the rest of this subsection, we will consider 
the momentum space to be a discrete set and consequently, the Dirac delta functions will be replaced by the Kronecker 
delta functions. As we shall see, the coefficients Cq „^ which determine the vacuum state of the full non-linear theory, 
are not all non- vanishing, by requiring that the exact vacuum state |0,i) satisfies the Schrodinger equation for the 
full Hamiltonian i.e. 

i^^\0;t) ^ (Hg + XHp)\0;t) (125) 

Using Eg. 1124(1 we get 

00 00 

E E E ^X'cg^.Jm,n2,..;t)G^J2 E E X'+'C^oL..Hp\ni,n2,..;t)G (126) 

/=0 ni,r!2.... fei,fe2,... ;=0 ni.n2,--- ki.kl,--- 

Comparing the coefficients of A' on both sides, we get G'^)i^^n2,..,nj,.. — <5o,ni<5o,n2----'^o,nj Then Ea. (|126|l reduces to 

00 00 

E E E *^'C'o!l....l^i>^2,..;t)G = E E E >^'ci;-^tMni,n2,.At)G (127) 

l—l ni,n2,... ki,k2,... l — l ni,n2,... ki,k2,... 

Equating the coefficients of the same power of A on both sides, we get the following equation for the expansion 
coefficients 



C'o!l,n2,...W = -« E E cl,,\]_^^,Jt)G{ni...n,..-t\Hp\m,,...rn,,..-t)G (128) 



The vacuum state correct to 0{X^) then becomes 

|0;i)[2] = |0;t)G + A ^ ^ Wlm, •■;i>G 

ni,r!2,---"j.--- ki,k2,...kj,... 

+ E E ciX...,n,Jt)\nu..,n,,.-t)G (129) 

ni,ri2,---"j.--- ki,k2,...kj,... 

From Ea. H123|l and Ea. (|128|) . the only non- vanishing contribution to Cq'^^^_^ ^. {t) comes from the following term 

l! j[dki][dk2\[dk'i][dki]5{ki + k2 + k^, + fc4)'^lfc^'^lfc2'i^'lfc3¥'lfc4aLfcjaLfc2aLfe3aLfc4, where the indices correspond to 
excitations in modes ki,k2,k3. The most general form of the equation for the first order coefficient follows from 
Ea. (|128l) . It can be written as 

where {n^j} = , n^.^, n^jg , nfc,^ and and Nf are numerical coefficients which depend on the relationship between 
the modes ki, k2, k^, k^. Nc correspond to the factor determined by the maximum power of a creation operator, 
whereas Nf is a combinatorial factor. It is important to note that distinct first-order coefficients (with distinct values 
of Tifcj^ , nk2 etc.) arise when constraints are imposed on the momentum modes ki,k2,k3, k^ which account for different 
types of scattering processes. In the absence of any constraints, Nc = Nf = 1 and Ea. p3U|l gives 

C'o;ll,,lfc2,lfc3,lfc4 = -^_^k,'Pk2fkVkJki+k2+ks+k^.,0 (131) 
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When the four modes appearing as indices in the creation operators are paired resulting in two distinct creation 
operators (i.e. when fcj = k'l and k'l = k'^), Nc = (-\/2)^ and Nf = {Ac^l'^) x 2c2 (since there are 3 ways of forming 
2 distinct pairs of creation operators); Ea. H130|l becomes 

C^ik,,2,, - -|(V2)'(¥'K)'4,+fe.,o (132) 

For a single pairing (i.e. when for example = k'l but k'l ^ k'^), Nc = \f2 and N f = (since there are ways 
of forming a single pair out of 4 creation operators) and Eq. (|130|) becomes 

From the form of the string of creation operators appearing in Ea. p3()|l . it is clear that the number of different 
coefficients correspond to the different ways of partitioning of four particle excitations in terms of excitations of lower 
order. 

In order to derive non-Gaussian corrections to the 2-point functions, it is necessary to obtain the corrected vacuum 
state, at least upto Oi}?). For this purpose, we need to determine the second order coefhcients. Although, there are 
many distinct second order coefficients, we will see that only a few of them will contribute to the 2-point functions, 
thereby considerably simplifying the calculations. Nevertheless, for the sake of completeness, we explicitly outline 
below, the general method for obtaining all possible second order coefficients. The equation for the second order 
coefficient can be written in the most general form 

It is clear from the string of annihilation and creation operators appearing in each term in Eq. (|123l) that the second 
order coefficients can be divided into five classes corresponding to zero, two, four, six and eight particle excitations 
respectively. The zero-particle excitation case leads to just one second order coefficient given by the equation 

/^(2) _ « ^(1) r 

(ni, ..,nj, ^ 4i'+fc^'+fc^'+fc^^o<<2fei'<<2fe^,'<<2^'<y2fei'afc;'afc^^^ (135) 

h'l ,/C2 ,■ ■ 

By making use of the commutation relations between the annihilation and creation operators to normal order the 
above string of operators, it can be easily shown that the above equation reduces to 

^ofo = -^C'^^U^4fc2.1k34fc4^'=l+'=2+fc34-fc4,0V'fcl</'fc2<<5fe3</'fc4 (136) 

The second order coefficients corresponding to the two particle excitation case, which is contributed by the second 
term in Ea. H123| ) can be written in general as 

(m, ^ ^-k'i+wi+w-+ki,Qh'^,k'ih-^.k'ih'^M'l'Pl^^^ (137) 

k'^ ,k^,.. 

which leads to the equations for the following second-order coefficient 

C'oafci,U5 = -^^^C^^li,^M^M^M^^^i+k.2+k-i+ki.n^k^+k^+ki~k^^ (138) 

The third term in Eq. H123|) leads to the following general equation for the second-order coefficients associated with 
the four particle excitation case 

(ni, ..,nj, ..■,t\ Y ^-k-~k'^+k'^+k'^.oh'^,k!jSk'^^k'^ipl,,(pl,_,(pk'^^^^^ (139) 
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which leads to the dynamical equations for the following four second-order coefficients 



(2) 



.(2) 



= -3iC,^;^i^^^ I^^(5fci+fc2+fc3 + fe4,043 + fe4-fe5-fe6,0<y5fc5'^fcg<,5fc3^ 

.(1) 



^0%k, ,lfc2 - -12i(\/6)Co^^i^^ ,lk^,lk3.1ki ^-2fel+fe3 + fe4,04i+fc2 + fc3+fc4,0'/'fei 'Pfc3<y^fc4 



(140) 



Similarly, the second order coefficients corresponding to the six and eight particle excitations contributed by the 
last two terms in Ea. H123|l respectively can be easily obtained. The general structure of the second order coefficients 
corresponding to the six particle excitation case is given by 



/v(2) _ _2i ^ (1) 



1 "^2 3 4 



1 "^^2"^"^^ ~^^4 



and that associated with the eight particle excitation case is given by 



-A(2) _ _ « ^(1) c 



1 "-2 3 4 



However, none of those coefficients contribute to the expressions for the two-point functions correct to 0{X^). Hence, 
we do not list them explicitly. 

The normalization of the vacuum state correct to second order in A is easily found in terms of the first and second 
order coefficients 



[2](0;i|0;i)[2] = l + A^ 



0;2fc^ ,2^2 0;2fc^ ,2^2 / ^ 0;2fc^ ,1^2 ,lfc3 0;2fc^ ,1^2 ,1^3 

fei,fc2 ki,k2,k3 



ki,k2.k3,k^ 



(143) 



We can now obtain expressions for the equal-time, 2-point functions by taking the expectation values with respect 
to this improved, albeit perturbative, vacuum state, with the understanding that 



(144) 



where the subscript indices a, 6, • • • = 1, 2 and $q is to be understood as the field operator ^{k) (a — 1) or the conjugate 
momentum operator H(fc)(a = 2). Moreover we will be working with the symmetric theory which implies that = 0- 
It is important to note that in the LvN formalism, the field and momentum operators are time-independent. The time 
dependence of the 2-point functions arise through the time-dependence of the auxiliary field variable ip. It is then 
easy to show that the equal-time, 2-point functions are given up to and including C'(A^) by the equations 



G?^\k;t) 



1 , \2/v^^(l)* ^(1) , ^(1)* ^(1) 

+ 1 '-'0;2fc^ ,2fc2 '-'0;2fcj ,2fc2 Z^ '-'0;2fcj ,lfc2 '1^-3 0^2^^ ,1^2 ^1 

fel,fc2 ki,k2,k3 



Er^W* /^(i) 4_ cr^(2)* ^(2)^ 

^0;lfci,U-2,lfc3,lfc4"-'0;lfci,lfc2,lfc3,lfc4 + l'-'0;0 + ^0;0 J 

fcl,fc2,fc3,fc4 

4 E (^0;2^1.2., G^1..2., V-k^-k + C«:.2., ^l^k 
ki ^ 



0/^(2) 



0(A3) 



(145) 
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G^?(fc;i) 



1 I \2/v^^(i)* ^(1) , ^(1)* 

E^(l)* /^(l) 4. Cr(2)* I /^(2)n, 



+ 4 E ^o;i:.,2., c«_^^2., ^-.^-^ + c^;i:2., c^^^^., ^i^^ 



+ A2 



ori(2) 



(146) 



G22(fc;0 = 



1 I \2/v^^(l)* ^(1) , ^(1) 

+ 1 '-'0;2fcj ,2^,2 0'2<«i .2fc2 ^ "-'0;2t, ^ 4fc2 > l^s 0;2fc ^ 4^2 > If^a 

fci,fc2 ki,k2.ks 



Z^ "^0;lfci,lfc2 4fc3,U4^0;lfci,lfc2.1'=3'lfc4^ ' 

fei ,k2,k3,ki 



(2)* , ^(2)- 
0:0 "t''-"0:0- 



4E(<21.2.,G^;2L.2.,^1;.-^-^ + ^^^1*.2.,^^^1,2.,^I-^.- 
fcl ^ 



+ A 



ori(2) 

^'-'0;l-k,lk'^-kVk 



(147) 



It is clearly evident from this set of equations that the first term corresponds to the Hartree approximation result 
and corrections to the Hartree approximation appear only at ©(A^). The auxiliary field mode variable (pk is given 
by the solution of Ea. (|101|) with A = 0. The effect of scattering enters through the presence of the first and second 
order coefficients which first appear at 0{}?). For example, 4-particle scattering processes with overall momentum 
conservation are encoded in terms involving the first order coefficients cl^^ i i i • The LvN formalism thus 
provides a systematic perturbative method for computing non-Gaussian (beyond Hartree) corrections to the 2-point 
function and higher n-point functions. Systematic corrections in powers of A > 3 can also be derived in a similar 
manner after including terms of C(A") (n > 3) in the expression for the improved vacuum state. These expressions for 
the 2-point functions particularly suited to obtain non-Gaussian corrections to the domain size in systems undergoing 
a quenched second-order phase transition. The 4-point function and higher n-point functions in momentum space can 
be also be computed in a similar manner, in terms of the auxiliary field modes and the expansion coefficients up to 
0{\^). The A-independent terms in the 4-point functions correspond to the factorization of the 4-point function in 
terms of products of 2-point functions and thereby represent the Hartree approximation result. 



B. Heisenberg Formalism : Evolution beyond the leading order 



In this subsection, we obtain the nonequilibrium evolution equations for the connected correlators, by taking the 
vacuum expectation value of the Heisenberg equations of motion for various combinations of products of field and its 
conjugate momentum operators, after making use of the cluster expansion to express the ordinary n-point functions 
as sum of products of connected n-point functions of lower order. This technique was used to obtain the effective 
potential and investigate phase transitions in spontaneously broken (^''-theory in 1 -I- 1 and 2-1-1 dimensions [s^l ■ We 
note from Sec. II that the expectation value of a functional of Heisenberg-operators with respect to the vacuum state 
becomes 

(0|/'H(i)|0) = (0,t|/'s|0,i), (148) 

where FY{{t) is any general function(al) of operators in the Heisenberg picture and i^s is the corresponding time- 
independent operator in the Schrodinger picture. 

From now on, to take into account NLO effects and compare with other methods, we adopt an approach different 
from the LvN formalism. We take the expectation value of the Heisenberg equation for any functional of (t^nit) and 
Trt{{t) as follows 

(0|z|/'H(0|0) = (0|[FH(t),^H(t)]|0) (149) 

where Hii{t) and Fait) are time-dependent operators and |0) is the time-independent vacuum state. Even though the 
explicit form of Huit) and FYi{t) in terms of the Schrodinger operators requires knowledge of the unitary evolution 
operator, we will not be required to know their explicit forms for the purpose of the calculation below. 
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In general, the evolution equations for the ordinary n-point correlators which are defined as 

5'i"l,;,(xir • • ,x„;t) = (0a(xi,i)---(^6(x„,t)) (150) 
= (-FH;a,---,b(xi,- • • ,x„;i)) 

are obtained by taking the vacuum expectation value of the Heisenberg equation (after dropping the subscript 'H') 



d{F,,..,t) I p 



,x„;i),i?]; 



(151) 



where the subscript indices a,b, - ■■ = 1,2 for our model Hamiltonian and (f>a is to be understood as the fluctuation 
field operator (j}f{'K,t) (a — 1) or the conjugate momentum operator 7r/(x, i) (a — 2) in the Heisenberg picture. The 
subscript / has been removed for notational convenience. The fluctuation field and its conjugate momenta satisfy the 
usual commutation relations 



[<^Q(x,t),(^6(y,i)] 



i6ab6^ (x 



y)- 



(152) 



For the unbroken symmetry case, (jjc = 0, only correlation functions of even order are non- vanishing. (For a theory 
exhibiting spontaneous symmetry breaking, a non-vanishing vacuum expectation value will induce a cubic interaction 
term because of which even odd n-point functions become non-trivial). Using the cluster expansion j37j | to express 
ordinary n-point functions in terms of the equal-time, connected n-point correlators, it is possible to obtain a set of 
evolution equations for the connected correlators. Because of the presence of the quartic coupling, the equations for the 
2-point functions would depend upon the connected 4-point functions; the equations for the 4-point functions would 
depend upon the 6-point functions and so on, thereby yielding an infinite hierarchy of evolutions equations for the 
n-point correlators. To go beyond the Hartree approximation requires including the cluster expansion of the 6-point 
function in terms of products of connected n-point functions of lower order and ignoring the effect of the connected 6- 
point function in the expansion. This amounts to taking into account direct scattering effects and provides a systematic 
way of going beyond the leading-order mean field expansion. This then yields a set of closed equations which have 
been truncated at the 4-point level. The leading order (Hartree approximation) result is easily recovered by expressing 
the 4-point function as products of connected 2-point functions and ignoring the term involving the connected 4-point 
function in the cluster expansion. This yields a closed set of equations for the equal-time, connected 2-point functions. 
In this scheme, incorporating NNLO effects would then amount to truncating the hierarchy of evolution equations 
at the 6-point level. Appendix A lists the cluster expansion of n-point functions (n < 4) and the set of evolution 
equations for the equal-time, connected n-point functions {n < 4) in configuration space. It is often more convenient 
to work with the Fourier transforms of the evolution equations. For that purpose, we define the Fourier transforms 
of the 2 and 4-point functions as 
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51^(^1,2:2) = y Mfci][dfc2]Gif)(fci,fc2)e*('=-^^+'=^-"^) 
alTJ {xi,x2,X3,X4) ^ I [dfci] [dk2] [dfca] [dk^] Gil2 (^1,^2,^3, fc4)e*(^ -^'^ 



For translationally invariant theories, 

Gif^(fci,fe) - G,('^^ 



t (fcl)(5(fci +fc2) 



Gil2 (fci , fc2 , A:3 , ki) = Gil2 (fci , ^2 , A:3)<5(fci + fc2 + fca + fc4) 



(154) 



which is just an indication of the conservation of momentum at the vertices. The Fourier transform of Eq. (Al) and 
Eqs. (A3)-(A7) yields the following evolution equations for G\^^ (fc): 



Glf^(fci) 
G^f^(fci) 



G^f^(fci) 



Gif^(fci) + Gri"^(fci) 



(2C), 



G^f^(fci) 



^'(fcl) 



[dk']G^^^\k' 



Gif^(fci) 



[dfc'][dfc"]GitS(fc',fc",-fci), 



+ 9 / Mfc']Glf)(fc') 



Gif^(fci) + G^f^(fci) 



j [dk'\[dk"] (^G(tS(fc', k", fci) + G(t2(fc', k", fci) 



(155) 
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By neglecting terms involving the connected 4-point functions in the above equations and setting (j)c = (since we 
are discussing the symmetric theory), one recovers the leading order (Hartree approximation) evolution equations, 
Ea ll()7l obtained in the last section. This also establishes the equivalence of the method used in this section to obtain 
NLO equations and the LvN method used in the previous sections. The evolution equations for the 4-point functions 

G[f,l{ki,k2, h) = Gif,l{ki,k2, fcs) + G^^2n{ki,k2, fcs) + G^t^^ (fci, fe, k^) + G[f,Uki,k2, fcs), (156) 
G'2m(fci,fe,fc3) 



mu me IjVin meLnou useu in tne previous seciions. ±ne evoiuiioi 
obtained by taking the Fourier transform of Eqs. (A3)-(A7) are 

Gif,l{ki,k2, fcs) + G^^2n{ki,k2, fcs) + G^t^^ (fci, fe, k^) + G[f,Uki,k2, fcs) 
G'22n(fci, fe, h) + G^t?i (^1. ^2, fcs) + G^ti2(A;i, fc2, k^) 
-uj^ki + k2 + fc3)Gltn'(fci, fc2, fcs) - AG^f )(fci)G(f' (fc2)Glf Hfca) 



[dk'] G 



;(f )(fc')Gitn'(fci, fc2, k,) + Glf (fcOG^tn^fc', fe, fcg) 



+G(f ) (fc2)G(tn^ (fc', k, , k,) + G^f ^ (fc3)G(tn^ {k', k, , k2) 



(157) 

G^tS(fci, fe, fcs) = G^t2(fci, fe, fca) + Gif.lik,, k2, fcg) - ^'(fci + k2 + k3)G[f,l{k,, k2, fca) - c^'(fci)G^tn^(fci, fe, fcg 
~x( G^f ^(fci)G(f ^(fc2)G^f' (fca) + G^f ^(fci)Glf ^(fc2)Glf' (fcg)' 



:dfc'] ^G^f )(fc')G(tn^(fci, fc2, A:3) + G[f\k,)G[f;i{k' , k2,k,) 



A 
'2 

HG^f ^ {k2)G[f,l (fc', fci , fca) + G^f ^ (fc3)Glt?i^ (fc', fci , fc2 
H J [dk'] (^G^f )(fc')G^tS(fci, fc2, ^3) + G^f ^(fci)G(tn^(fc', fe, fcs) 
l-G^f )(fc2)G(tS(fc', A:3) + G^^,''\k,)Gif,l{k' , k,, k2 



(158) 



G2t2l(fcl,fe,fc3) 



G't2^l{ki,k2, ks) - uj^{ki + k2 + k'i)G'f^l{ki,k2, ks) 
-u;2(fci)G(tS(fci, k2, fcg) - u;'(fc2)G(tn^(fci, fe, fcg) 

- A (^G^f ^ (fci)Gif' (fc2)G(f ) (fcg) + G^f ) (fci)G(f ) {k2)G[f^ {k,) 
+Gif\k,)Gif\k2}G[f\k,)' 



"21 

J [dk'] (G[f\k')G[f,l{k,,k2, ks) + G[f\k,)G[f,l{k', k2,ks) 
+G(f' (fc2)G(t?^ (fc', k, , A:3) + G^f ) (fc3)Glt?2^ (fc', fci , fc2)) 
+ 1 [dk'] (^G^f )(fc')G^t?^ (fci, fc2, A:3) + Gg^\k,)G[f,l{k', k2, fcs) 
+G(f' (fc2)G(tS ik', k, , A:3) + G'f ) (A:3)G(tS (^'^ , ^2)) 
+ J [dk'] (^G^f \fc')G^tn'(fci, fc2, ks) + G^f )(A;i)G(tn'(fc', fe, fcs) 
+G(f' (fc2)G(tn' (A;', A:i , ^^3) + G^f ^ (A;3)G(t?i' (fc', fci , fc2) 



(159) 
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+u;^ik2)Gg^,lik^,k2,k,) + c.2(A;3)G(tS(fci, fc2, k,)^ 

If )(fci)G(f Hfc2)G(f )(fc3) + G^f )(fci)G(f )(fc2)G(f )(fc3) 
if )(fci)G(f )(fc2)Glf )(fc3) + G^f \fci)G(f ^(fc2)G(f )(fc3) 



-A G 



+G5, 

_A 
~2 

+G 



if Hfc2)G(f2^(fc', fci, fc3) + Gff '(fc3)G(t2(A:', fe)' 



[dk'] ( G 



hG 



hG. 



if )(fc')G(tS(fci, fe, fca) + G^f )(fci)G(t2(A;', fc2, ks) 
If ^ (fc2)G(tS (fc', fci , fcs) + Glf ) (fc3)G^f 2^ {k', k, , A;^)) 
/' [dk'] (G^^f\k')Gif,l{kuk2, fca) + G^f )(fci)G(f k,,k,) 
ff\k2)Gif,l{k', fci, fcg) + Gif )(fc3)G<f 1^ (fc', kuk^)^ 
f [dk'] (G'^;f\k')Gif,l{kuk2, k,) + Gif\k,)Gif,l{k', k,,k,) 
if )(fc2)G(f )(fc', fci, fcs) + G^f )(A3)G^f (fc', A;^)^ 



(160) 

This set of equations completely determine the evolution of the correlation functions beyond the leading order. 
Numerical integration of this set of equations with noncquilibrium initial conditions and tracking of the subsequent 
dynamical evolution of the connected, n-point functions would make it possible to ascertain whether the proposed 
truncation scheme is good enough to ensure late-time thermalization of the system. The inclusion of the 4-point 
correlators in the hierarchy amounts to considering terms of 0{hP) in the effective action. Since each connected n- 
point function G^"*^-* ~ 0{h'^~^), inclusion of higher order connected n-point correlators in the hierarchy of evolution 
equations allows for a systematic way of incorporating quantum effects. The issue of truncation error arising from 
truncating the hierarchy at the 4-point level needs to be addressed. This issue has been briefly discussed in Ref. 
where it was argued that truncation error would build up with time and eventually invalidate the working assumption 
of the formal hierarchy of the connected correlators given above. For the quantum mechanical anharmonic oscillator 
and the 0(N) vector model, the decoherence time-scale for breakdown of the hierarchy was qualitatively argued to 
scale with h~^^'^ and ^/N respectively. However, their analysis was based on a simple (0 -I- 1) dimensional quantum 
mechanical model and quantitative treatment of this aspect for a system with infinite degrees of freedom (field 
theory) will certainly be more complicated. Although, inclusion of higher order n-point functions in the hierarchy 
does not eliminate the instabilities, it does postopone the onset of these instabilities. This delay in the appearance 
of intsabilities would make it possibile to gain valuable insights into the non-equilibrium dynamics and approach to 
thermalization, in a temporal domain which lies beyond the domain of validity of the Hartree approximation. A 
possible way to avoid the instabilities arising due to the truncation of the hierarchy might involve adapting the Direct 
Interaction Approximation(DIA) developed by Kraichnan in the context of turbulent fluid dynamics. We are 
currently exploring this possibility. 

VII. COMPARISON WITH ALTERNATIVE METHODS 

Several approaches have been developed to understand the dynamics of fields in noncquilibrium field theory. It 
would be useful to make a comparative study of the different approaches, specifying the relationship between the 
different methods and the relative merits and demerits of each approach. In this section, we will address this issue by 
comparing our approach with other methods which have been used to obtain evolution equations for the equal-time 
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correlation functions. We will compare our canonical approach with two specific approaches [241 136| that have been 
recently discussed in the literature. 

Wetterich's method 36] is based on deriving an evolution equation for the partition function (or the generating 
functional for IPI graphs). The time evolution of the n-point correlation functions (or n-point vertex functions) is 
completely determined by the time evolution of the generating functional and the microscopic dynamical equations 
of motion. We will show below that the evolution equations for the equal-time correlation functions obtained using 
Wetterich's method match exactly with those derived in this paper using a canonical approach. 

For a (/i*-theory without spontaneous symmetry breaking, the evolution equation for the generating functional can 
be cast in the form |3a 



dtZ[j{x),h{x);t] = {Cci + Cqm)Z[j{x),h{x);t] 



(161) 



where Cci and Cqm are the classical and quantum parts of the Liouvillian operator respectively and are given by 



dx 



5h{x) 



h{x){ V 



5j{x) 



X 6^ 



dj{x) 6 6j{xy 



A 
4! 



dx 



h^x) 



6j{x) 



(162) 



Here j{x) and h{x) are source terms for the field and conjugate momenta, respectively, and Z[j{x),h{x)]t] is the 
generating functional for n-point functions, Z = Tr[exp(J c?x(j(a;)0(a;) -f h{x)ir{x)y)p\, where p is the density matrix. 

It is fairly straightforward to obtain the corresponding evolution equations for the generating functional for con- 
nected graphs W[j{x),h{x);t] and the generating functional for IPI graphs T[ipc{x),T:c{x)]t]. {(pc{x) and tTc{x) are 
vacuum expectation values of the field and conjugate momentum in the presence of sources j{x) and h{x), respectively). 
Z[j{x),h{x);t],W[j{x),h{x);t] and r[(fc{x),Trc{x);t] are related by the equation W[j{x), h{x);t] = In Z[j{x), h{x);t] 
and the Legendre transform T[ipc{x),iTc{x);t] = —hiZ[j(x),h(x);t] + J dx[j{x)Lpc{x) + h{x)TTc{x)]. The evolution 
equation for W[j{x), h{x);t\ which completely determines the evolution of the connected n-point functions then be- 
comes 



dtW[j{x),h{x);t] J dx 
X 

+ 4! 



jix) 



dx 



Sj{x)^ 



, 5W 6^W 
Sj{x) 6j{xy 



5W 
Sj{x) 



3> n 



h^x) 



dW 



For a symmetric theory ipcix)\j=Q 



Sj{x) 
(. <5iv 



(163) 



•(2C)/ X 
5il (Xl,X2) 



.(2C)/ X 
921 [Xl,X2j 



.(2C), N 
522 {Xi,X2) 



5j{xi)5j{x2) 
6h{xi)Sj{x2) 



6h{xi)6h{x2) 



dtW[j,h;t]\j^h=a 
dtW[j,h;t]\j^h=o 



and 7r,(x)U=o = {^))\h=o = 0. Using Eq. we find that 

= .9l2 (a;i,2;2) +521 {Xi,X2), 

= 522 [Xi,X2j+yigi2 '[Xi,X2j-m [Xl,X2j 

A (2C)/ N (2C)/ N ^ (iC), X 

-o9ii [xi,X2)gii '(xi,xi) - -gl-^^{{xi,xi,xi,X2), 
2 6 



dtW[j,h;t]\j^h=o = (V?-m2)5(f)(xi,X2) + (V2-m2)g(f)(xi,X2) 



Xl {2C)f ^ (2C)/ N , (2C)/ N (2C)/ N 

2(511 (a;i,a;i)5i2 (a^i,a;2) +5ii (a:2,a;2)52i (2^1,2:2) 



5i 



1112 (2^1, a;i, Xl, 2:2) + g'2ii{(Xi,X2,X2, X2j 



(164) 



The above set of equations for the connected 2-point function, derived using Eq. H163|) is identical to Eq. (Al) which 
was derived using canonical methods and the cluster expansion for ordinary n-point functions. Similarly, one can 
show that the evolution for the connected 4-point function obtained using the above method is the same as the ones 
obtained in Appendix A. 

More recently, Ryzhov and Yaffe (RY) developed another method for obtaining the nonequilibrium, coupled set 
of evolution equations for the n-point correlators (24l | . Their method is based on an expansion of the coherent state 
expectation value of the products of operators in terms of subtracted n-point functions and the assumption that the 
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initial state is some coherent state. The central idea behind the RY method is the expansion of the Hamiltonian 
operator in terms of powers of the generators of the underlying coherence group, which is the Heisenberg group 
for the (j)^-theoTy. It consists of the operators e'-'^', e*^^), e^'^^ = $1 (fc), $2(fc), 1 satisfying the commutation relations 
[e^"'^-', e^^-*] = i/j^26^^-''5(fc — k') where /fg = — /li = 1 are the only non-vanishing structure constants. By taking 
the coherent- state expectation value of the Heisenberg equations of motion for appropriate products of equal-time 
coherence group generators, it is possible to obtain the evolution equations for the subtracted n-point functions (or 
connected functions) in the same manner as described in Sec. V. In Appendix B, the evolution equations for the 
two-point function are derived in the Hartree approximation, using the RY method and are shown to be equivalent 
to those obtained using the LvN approach. 



VIII. CONCLUSION 



In this paper we have used the LvN formalism, a canonical method, to study the noncquilibrium evolution for the 
equal-time, connected n-point functions for a symmetric (/)^-field theory. The usefulness and simplicity of the LvN 
method in obtaining perturbative, non-Gaussian corrections to the n-point correlation functions were first illustrated 
using the quantum-mechanical anharmonic oscillator model. The formalism was then applied to a (/)^-field theory 
and used to obtain the evolution equations for the connected 2-point function in the Hartree approximation, for 
noncquilibrium evolution as well as for the thermal equilibrium case. We also used the LvN formalism to go beyond 
the Gaussian approximation and obtain expressions for the 2-point functions correct to 0{X^) after calculating the 
improved vacuum state. Expressions for 4-point and higher n-point functions can be similarly obtained after some 
straightforward but tedious algebra. The corrections to the Gaussian approximation were found to appear first at 
©(A^). The noncquilibrium evolution equations beyond the leading order approximation were then obtained by taking 
the vacuum expectation value of the Heisenberg equations of motion for appropriate products of field operators. This 
provides an alternative and non-perturbative approach to investigate systematically the noncquilibrium evolution of 
quantum fields and yielded an infinite hierarchy of coupled equations for the connected n-point correlators which were 
truncated at the 4-point level. This involved ignoring the contribution of connected n-point functions (n > 6) in 
the evolution and resulted in a closed set of equations involving connected 2-point and 4-point functions. Since we 
restricted our investigation to a symmetric (/i^-theory with vanishing vacuum expectation value, all connected odd-n- 
point functions vanish and have no effect on the evolution. We have also established a connection between the canonical 
approach used in this paper and other methods |24ll36| developed in the literature for deriving noncquilibrium evolution 
equations for equal-time correlation functions. It would be straightforward to generalize this technique for studying 
the noncquilibrium evolution of spontaneously broken field theories. Symmetry breaking would result in the generation 
of linear and cubic terms in the potential, as a consequence of which, even connected odd-n-point functions would 
contribute to the dynamical evolution. 

The inadequacy of the Hartree approximation in studying the approach to thermalization has been extensively 
discussed in the literature. The Hartree approximation neglects scattering effects and therefore cannot account for 
long-time thermalization of the system. An elegant interpretation of this aspect has been put forward by Wetterich [3^ 
and is based on the realization that the Hartree solution corresponds to a fixed point of the theory (i.e. configurations 
for which dtZ = 0). The presence of an infinite number of conserved quantities (fixed points) prevents the system 
from escaping from these noncquilibrium fixed points and approaching the thermal equilibrium fixed point unless the 
initial values of all these conserved quantities coincide exactly with the ones corresponding to a thermal distribution. 
To take into account the effect of scattering which would ultimately lead to thermalization of the system, requires 
going beyond the Hartree approximation. An approach based on the loop expansion of the 2PI effective action [l^ 
in powers of h (or for an 0(N) symmetric field theory) holds more promise in this context. In this paper we 
have discussed two different canonical approaches (perturbative and non-perturbative) of going beyond the Hartree 
approximation. The perturbative LvN approach developed in this paper provides an elegant method for obtaining 
corrections to the n-point functions in terms of powers of the coupling constant A. The effect of scattering which first 
appears at O(A^) is manifest through the first and second order coefficients whose dynamical equations were derived. 
This method provides a straightforward method for obtaining non-Gaussian corrections to the domain size in systems 
undergoing a quenched second order phase transition. By making use of the exponentially growing solution of the 
soft modes of the theory, it is possible to obtain expressions for the first and second order coefficients, which in turn 
lead to expressions of the 2-point functions to various orders in A. The Fourier transform of the 2-point functions 
would then yield the domain size which includes corrections due to non-Gaussian effects Moreover, the LvN 

method provides an analytical, albeit perturbative, method for studying the role of interactions and effect of scattering 
on thermalization of the system. This is possible by obtaining perturbative corrections to the 2-point functions to 
various powers of A and comparing the resulting expression with the thermal 2-point correlation functions. We 
have also discussed in detail an alternative, non-perturbative canonical approach based on the Heisenberg formalism 
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to study the nonequihbrium evolution of the (j) field theory. This non-perturbative Heisenberg formalism requires 
incorporating the effect of connected 4-point function in the nonequihbrium evolution equations. The effi cacy of 
truncating the hierarchy at the 4-point level in ensuring thermalization needs to be addressed. Recent work |23| for 
classical field theories has been inconclusive because of the prohibitively long time required for equilibration and also 
due to the fact that the numerical evolution becomes unstable long before the system thermalizes. We believe that 
the canonical approaches developed in this paper provide an useful alternative method for studying nonequihbrium 
dynamics in field theory. 
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APPENDIX A 

In this appendix we write down the set of evolution equations for the connected n-point functions in configuration 
space by making use of the cluster expansion for the 4-point and 6-point functions. The cluster expansion of n-point 
functions can be formally derived from the generating functional of the interacting field theory |37j . The difference 
between the ordinary and connected n-point correlators first appears at the 4-point level. To simplify the notation, we 
represent the (j) and tt field operators by the numbers 1 and 2 respectively. So the correlation functions g^^ {xi,X2), 

9^^^\xi, X2) etc. are symbolically represented as gj^^'^'' (1, 2), g^^'^\l,2), respectively, where the numbers in the 
brackets refer to the subscript indices of the spatial coordinates. 

The equations for the 2-point functions in configuration space are easily obtained by using Eq. (|151|l and the 
appropriate cluster expansion of the 4-point functions. They are 

Sjf) (1,2) = g[f\l,2)+g[f\l,2), 

5^^(1,2) . (V? -m^)5|f ^(1,2) + iVl-m^)gir\l,2) - ^ (5^^!, 1)5^^1, 2) + ^^^(2, 2)5^^(1, 2)) 

(s^tS (1, 1, 1, 2) + g^tn' (1, 2, 2, 2)) , (Al) 

where the last two equations in the above set have been obtained by making use of the following cluster expansion of 
the 4-point functions 

.g(tli(l,l,l,2) = 5lt2(l,l,l,2) + 3.9lf)(l,l)glf)(l,2), 
.gitl2(l,l,l,2) = 5itS(l,l,l,2) + 3.9if)(l,l)gif\l,2), 

5^tli(l,2,2,2) = .9^tnHl,l,l,2) + 3.9if)(2,2)g(f)(l,2). (A2) 

We have used lower case letters to represent the correlation functions in configuration space and upper case letters 
to represent their Fourier transformed counterpart in momentum space. The n-point connected correlation functions 
appearing above and in the subsequent discussion are all normal ordered vacuum expectation values of products of 
operators. 

The evolution equations for the connected 4-point functions in configuration space can be similarly derived by 
making use of the appropriate cluster expansion of the 6-point functions. The five independent equations for the 
4-point function are 

g^tn^l, 2, 3, 4) = s'tn (1, 2, 3, 4) + gil'^ln, 2, 3, 4) -t- ^(^^^^(l, 2, 3, 4) + ^(^^(l- 2, 3, 4), (A3) 
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5^tS(l, 2, 3, 4) = g^^lil, 2, 3, 4) + g'iglil, 2, 3, 4) + ^^tSCl, 2, 3, 4) + (V? - m')5itS(l, 2, 3, 4) 
-A 2),f,^)(l, 3),f,^)(l, 4) - I 

+5if ^(1, 3)5it2(l, 1, 2, 4) + ffif 4)5it2(l, 1, 2, 3)) , (A4) 

ij^^l (1, 2, 3, 4) = g'i^lil, 2, 3, 4) + g'^^l (1, 2, 3, 4) + (V? - m^)g^^l (1, 2, 3, 4) + (V^ - m')gt^l (1, 2, 3, 4) 
-A f fflf ) (1, 2)5if ^ (1, 3)5lf ' (1,4)+ gf^^ (1, 2)5if ' (2, 3)5if ' (2, 4)' 



-y (sif ^(1, 1)5^2^(1, 2, 3, 4) + gg^\l, 2)i?(tS(l, 1, 3, 4) 
+ff(f ^(1, 3)ffit2i^(l, 1, 2, 4) + ffif ^(1, 4)5it2i^(l, 1, 2, 3) 
- y (ff^^ (1, 2)fl^tS (2, 2,3,4) + ^^f ^ (2, 2)g^^l (1,2,3,4) 

+5(f )(2, 3)3(^^(1, 2, 2, 4) + ^(2, 4)5^tn^(l, 2, 2, 3)) . (A5) 

Since the correlation functions are all normal ordered, the above equation is symmetric under the interchange of 
indices x\ ■«-> X2 and X3 ^ X4 

#22(1, 2, 3, 4) = 5(^22 (1, 2, 3, 4) + (V? - m^)g^l{l, 2, 3, 4) + (V^ - m^)g^^,^l{l, 2, 3, 4) 
+(Vl - m^)^^^^ (1,2, 3, 4) - A (^g^^^\l,2)g^^^\l,Z)g^^,^\lA) 

+g(f ) (1, 2)g(f ^ (2, 3)5if ^ (2, 4) + gg^\l, 2,)gg^^ (2, 3)5if ^ (3, 4) 

- Y (^if ^(1, l)5i22i'(l, 2, 3, 4) + ff(f ^(1, 2)5it2(l, 1, 3, 4) 

+5^f '(1, 3)5ltS(l, 1, 2, 4) + ffif^l, 4)5^12(1, 1, 2, 3) 

^ ^ff^f \l, 2)5^t2(l, 2, 3, 4) + ff(f )(1, 2)5it2(2, 2, 3, 4) 



2 

+.g(f Hi, 3)fl^t9 (1, 2,2,4)+ 5(f) (2, 4)5^12 (1,2,2,3) 
- y (sif ^ (3, 3)5^^22^ (1,2,3,4)+ 5^1^) (1, 3)5^tfi' (2,3,3,4) 

+5^1''^ (1, 3)#i2 (1, 3, 3, 4) + (f ) (3, 4)5^^22 (1, 2, 3, 3) V (A6) 
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The above equation is symmetric under the exchange of coordinate indices xi X2 ^ x^. 

g^l^^il, 2, 3, 4) = (V? - m')g[l^i{l, 2, 3, 4) + (V^ - m^)gil^i{l, 2, 3, 4) 
+ (V2 - m')gif,^{l, 2, 3, 4) + (V^ - m')gg'^l{l, 2, 3, 4) 

-A (^ffif 2)g(f 3)glf 4) + gif\l, 2)g(f ^(2, 3)g[f\2, 4) 
+5^f Hi, 3)5^1''^ (2, 3)gf,^H3, 4) + g(f ) (1, 4)g(f ) (2, 4)g(^^^^ 



A 



.gif ^1, l)5i2S(l, 2, 3, 4) + 5if ^1, 2)5ilS(l, 1, 3, 4) 



+g(f 3)g(t2(l, 1, 2, 4) + g^f ^1, 4)ff(tS(l, 1, 2, 3) 
^ f.g(f ^1, 2)5it2(2, 2, 3, 4) + ^if ^2, 2)5^t2(l, 2, 3, 4) 



2 

+g(f '(2, 3)g^tS(l, 2, 2, 4) + g^f ^(2, 4)g(12(l, 2, 2, 3) 

-| (s^^l, 3).9^t?2'(2, 3, 3, 4) + .g(f )(2, 3)gif,^il, 3, 3, 4) 

+g(f )(3, 3)gif;^{l, 2, 3, 4) + g[f'\3, 4)gif;i{l, 2, 3, 3) 
A 



5^^!, 4).9^tn'(2, 3, 4, 4) + gif\2, 4)gif,l{l, 3, 4, 4) 
+5ri"H3, 4)g(^2S(l, 2, 4, 4) + g^f ^(4, 4)ff(^2S(l, 2, 3, 4)) . (A7) 



APPENDIX B 

We now show the equivalence between the LvN method and the RY method by obtaining the evolution equations for 
the Fourier transform of the subtracted 2-point correlations using the RY method in the Hartree approximation. The 
expectation value of the Hamiltonian, (|79|) . can be expressed in terms of the Fourier transforms, $oi(^7 0, = 1, 2), of 
the coherent state expectation value of the fluctuation field (j)o{x, t) = (0/)cs and that of conjugate momentum operator 
7ro(x, t) = (7r/)cs, respectively. In order to apply the RY method, it is necessary to isolate from the expectation value of 
the Hamiltonian, the parts which depend only on the Fourier transforms ^oi{k,t) = ^o{k,t) and $02 (fc,^) = no(fc,t) 
of the coherent state expectation value of the field and its conjugate momentum operator. Using the relations 
(7r^(a:))cs — 7''o(^,*) + ff22(0,t) and (0j(a;))cs — (j)Q{x,t) + gii(0,i), and the Hartree factorization of the cubic and 
quartic terms specified in Eq. (|94|l . the coherent state expectation value of the full Hamiltonian can be expressed as 

(i?(t))cs = J[dk]H{<i>o{k, t), Ho(fc, <)) + •••, (Bl) 

where the dots correspond to terms which are independent of <i>oi(fc, t), (i — 1,2), and 



1, 



2 



A .0 _ A 



ff(<i>o(fc,i),Ho(fc,t)) = -H^(fc,t)+ u;i + -^ii{t) + -gn{0,t) (B2) 



2""' ' 2- 



By noting that the underlying coherence group in this case is the Heisenberg group whose generators e^^-*, e^^\ e'^' 
= 7i'fe,<^fe,l satisfy the commutation relations [e'"'^' , e^^^] = ifi2G'^'^^5{k — k') where fi2 = — /li = 1 a-^e the only 
nonvanishing structure constants, and by making use of the general form of the NLO evolution equations for the 
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1-point and 2-point functions, we get the following self-consistent set of equations to 0{{X )^) 



dt 

dUQ{k,t) 
Jt 

dGii(fc,<) 
Jt 

dG22ik,t) 
dt 

dGi2{k,t) 
dt 

where we have used H'^^^'> = i- 



no(fc,i), 



A 



Gu{k,t) 



A 
2' 



2 

G2iik,t), 

'liit) + ^511(0, t)] (Gi2(fc,t) + G2i{k,t)) + 0((A 



A 



A 



G22{k,t) - [oji + -c}>ii{t) + -5n(0,t)jGn(fc,t) + 0((A )^), 
5[^' + l'^c/W + i5ii(0,t)] andil(22) ^ ^. 



(B3) 



l^.ak(k) = 5[^' + l'^c/W + i5ii(0,t)] andil(22) ^ ^ ^^r^tj^ = \. The set of equations 
(|B3|l are identical to the equations (|107(l obtained using the LvN formalism which clearly establishes the equivalence 
of the LvN and RY methods. The coherent state can be considered as the vacuum state of the theory with its 
expectation value providing the classical background. For a symmetric theory, the coherent state expectation value 
of the field and the conjugate momenta vanishes, and we recover the set of equations for the 2-point functions (in the 
Hartree limit) obtained in Sec. V. 
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